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We present a set of boundary conditions for solving the elliptic equations in the Initial Data 
Problem for space-times containing a black hole, together with a number of constraints to be satisfied 
by the otherwise freely specifiable standard parameters of the Conformal Thin Sandwich formulation. 
These conditions altogether are sufficient for the construction of a horizon that is instantaneously 
in equilibrium in the sense of the Isolated Horizons formalism. We then investigate the application 
of these conditions to the Initial Data Problem of binary black holes and discuss the relation of our 
analysis with other proposals that exist in the literature. 



PACS numbers: 04.20.-q, 04.25. Dm, 04.70.Bw, 97.80.-d 



I. INTRODUCTION 

The problem of determining appropriate initial data 
for binary black holes is of crucial importance in order to 
construct successful numerical simulations for these as- 
trophysical systems pj . Starting with Einstein field equa- 
tions, a specific strategy for this problem consists in solv- 
ing the relevant elliptic equations on an initial Cauchy 
surface where a sphere S has been excised |2| for each 
black hole (inner boundary) . The purpose of the present 
work is to present a set of inner boundary conditions in- 
spired by purely geometrical considerations, inasmuch as 
they are derived from the formalism of Isolated Horizons 
S [J El and which guarantee that the excised 
sphere is in fact a section of a quasi-equilibrium horizon. 

A pioneering work on this inner boundary problem was 
presented by Cook in Ref . , in the context of formu- 
lating a definite full prescription for the construction of 
initial data for binary black holes in quasi-circular orbits. 
The assumptions made in that analysis permit to deter- 
mine a proper set of conditions for a quasi-equilibrium 
black hole. However, at least intuitively, one would 
expect that the Isolated Horizons formalism, which is 
mainly a systematic characterization of the notion and 
properties of quasi-equilibrium horizons, could supply a 
more powerful and consistent framework for discussing 
the conditions in this black hole regime. Actually, the 
spirit in Ref. [jj closely resembles that encoded in the 
Isolated Horizons scheme, but does not fully capture it. 
Therefore, in this specific sense, the quasi-equilibrium 
horizon analysis may be refined. With this motivation, 
we will truly adopt here the Isolated Horizons formalism 
as the guideline of a geometrical analysis whose ultimate 
goal is the ab initio numerical construction of an isolated 
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horizon. This strategy provides us with a rigorous math- 
ematical and conceptual framework that systematizes the 
physical assumptions. 

For the sake of clarity, we have considered important 
to provide a relatively self-contained presentation, even 
at the cost of lengthening the article. The rest of the 
work is organized as follows. As in Refs. [HEl, we use 
a Conformal Thin Sandwich approach to set the Initial 
Data Problem; thus Sec. II briefly reviews the basics 
of this approach. Sec. Ill introduces the main ideas of 
the Isolated Horizons framework and underlines the im- 
portance of its hierarchical structure by first introducing 
Non-Expanding Horizons and then Weakly Isolated Hori- 
zons. Boundary conditions on the horizon are derived in 
Sec. IV. Sec. V discusses the relation of this approach 
with that of Ref. Q- Finally Sec. VI presents the con- 
clusions. 



II. CONFORMAL THIN SANDWICH 
APPROACH TO INITIAL DATA 

In this section we formulate the problem that will be 
analyzed in this work and introduce our notation. We 
will use Greek letters (/x, v, ...) for Lorentzian indices, in- 
termediate Latin letters ■■■) for spatial indices on a 
Cauchy slice, and Latin letters from the beginning of the 
alphabet (a, &,...) for coordinates on a two-dimensional 
sphere S 2 . 

Adopting a standard 3+1 decomposition for General 
Relativity (see e.g. Ref. ^H), the space-time M. with 
Lorentzian metric <? M „ is foliated by spacelike hypersur- 
faces E t parametrized by a scalar function t. The evolu- 
tion vector < M , normalized as i M V M i = 1, is decomposed in 
its normal and tangential parts by introducing the lapse 
function a and the shift vector /3 M 
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where n,, 



-aV^i is the unit timelike vector normal to 



E t and n"/3„ = 0. 
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Denoting by 7^ the induced metric on S t , the Lorentz- 
ian line element reads 

ds 2 = g^ v dx^dx v 

= -a 2 dt 2 + ^{dx 1 + (3 l dt){dxi + (3 j dt) . (2) 

The embedding of the hypersurfaces S t in the four- 
geometry is encoded in the extrinsic curvature: 



which can also be expressed as 

K ij = -±-(d t 'y ij -D i l3 j -D i p i ) 



(3) 



(4) 



where Di is the connection associated with 7^. 

Under the 3+1 decomposition, Einstein equations split 
in two sets: evolution and constraint equations. In vacuo, 
the case that we are interested in, the evolution equations 
are 

d t K i:j - C K i:j = a (R VJ ~ 2K lk K k 3 + KK i3 ) 

- DiDjd , (5) 

where K is the trace of Kij (K — j^Kij). On the other 
hand, the constraint equations (respectively, Hamilto- 
nian and momentum constraints) are expressed as 



KijK %: > = 



R + K 

Dj (K lj - f j K) = 



(6) 
(7) 



In brief, the Initial Data Problem consists in providing 
pairs (jijjK^) that satisfy the constraints © and Q on 
an initial Cauchy surface S t . 

The discussion of the notion of quasi-equilibrium de- 
mands a certain control on the time evolution of the rele- 
vant fields. The Conformal Thin Sandwich (CTS) intro- 
duced in Refs. |12l is particularly well suited, since 
it provides an approach to the Initial Data Problem that 
consistently incorporates (a part of) the time derivative 
of the metric, together with the lapse and the shift. 

The CTS approach starts by conformally decompos- 
ing the metric and the extrinsic curvature, the latter ex- 
pressed in terms of its trace K and a traceless part A %: > , 



lij = * lij i 
/\ ■' = # l A ij + 

3 



(8) 



In this expression the conformal factor ^> is given by 



* = 



(9) 



where 7 is the determinant of 7^ and / is the determinant 
of fij, an auxiliary time-independent metric, dtfij = 0, 
which captures the asymptotics of 7^ [T^j . 



Substituting the above decomposition of 7^ in relation 
I0J and taking the trace, we find: 1 

dt* = p i b i ^ + %{D i f3 i -aK). (10) 
6 

This expression will play an important role when setting 
the appropriate boundary conditions in Sec. IV. 

With Eq. (jHJ, the Hamiltonian constraint © can be 
written as an elliptic equation for the conformal factor: 



~ ~ . * - - / 1 - K 2 

DiD % ^> = —R — \& I -An A* 3 . . 

8 \8 13 12 J ' 



(11) 



whereas the momentum equation is expressed as an el- 
liptic equation for the shift 



1 ^ .-, 



DjD 3 ^ + -DWjfi + R i 3 l3 j - (Lpy j Dj ln(a*" 6 ) 
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Djdtf 3 +d t f J D j ]n(a^- 



(12) 



where Aij = "(ikljiA M , Di is the connection associated 
with and (Lfi) l i=D l l3 j + Z) J f3 l - §Z) fe /3 fc 7 y '. 

Only the conformal part 7^ of the metric 7^ encodes 
dynamical degrees of freedom. This suggests solving the 
trace of the evolution equations JSJ together with the 
constraints, defininghi this way an enlarged problem on 
the initial surface |l3j |. This additional equation turns 
out to be elliptic in its dependence on the lapse a, 



D i D' l a + 2Di In* 17a 
= * 4 la ( AhA* + 



P l DiK - d t K \. (13) 



K- 



An extra justification to add this equation is because 
it straightforwardly permits to impose the condition 
dtK = 0, a good ansatz for quasi-equilibrium. In this 
extended problem the constrained parameters are given 
by f3 l , a) and the free data on the initial Cauchy sur- 
face are (7^, 9t7 y , K, dtK), subject to the constraints 
det(7ij) = / and Jijdtj 11 = [in the strict Initial Data 
Problem ^ is a free parameter on the initial slice, but 
here it is constrained owing to Eq. I|13|) ]. Hence, the inner 
boundary problem presented in the Introduction reduces 
to the search for appropriate boundary conditions for <J/, 
/3 % , and a imposed on the horizon. 



III. ISOLATED HORIZONS FORMALISM 

In this section we will motivate the introduction of 
the notion of Isolated Horizon (IH) and summarize the 



1 In the following expressions, objects with an over-tilde are asso- 
ciated with the conformal metric ^ /g. They are consistent with 
the conformal rescalings in Ref. JTJ , rather than with those 
originally introduced in Refs. Il2tll I . 
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concepts and definitions that will be employed in Sec. 
IV. We will try to provide a presentation as accessible as 
possible for a broad community. For a detailed and more 
rigorous discussion of the IH formalism see, e.g., Ref. Q 
or Ref. 

The physical scenario that the IH construction at- 
tempts to describe is that of a dynamical space-time con- 
taining a black hole in equilibrium, in the sense that nei- 
ther matter nor radiation cross its horizon. This scenario 
applies as an approximation for each of the two black 
holes in a binary before their merger, provided that they 
are sufficiently separated, therefore justifying the rele- 
vance of the IH formalism for the Initial Data Problem 
of binary black holes. 

A very important feature of the IH formalism is its 
(quasi-)local character. In our context, the need of a 
(quasi-)local description is motivated first by the way 
in which numerical simulations are designed from a 3+1 
approach, in which we do not have a priori control on 
global space-time properties, and secondly by the desire 
of characterizing physical parameters of the black hole 
as well as the concept of equilibrium in a (quasi-)local 
manner. The notion of apparent horizon, with a local 
characterization as an outermost marginal trapped sur- 
face 2 in a three-slice, seems an adequate starting point. 
However, to include the concept of equilibrium we must 
somehow consider the evolution of this two-dimensional 
surface. In the (quasi-)equilibrium regime, the notion of 
the world-tube of an apparent horizon does in fact make 
sense (there are no jumps). Actually, an IH implements 
the idea that an apparent horizon associated with a black 
hole in equilibrium evolves smoothly into apparent hori- 
zons of the same area, in such a way that the generated 
world-tube is a null hypersurface. This null character 
encodes the key quasi-equilibrium ingredient, and is es- 
sentially linked to the idea of keeping constant the area 
of the apparent horizon. 

Inspired by these considerations, the definition of IH 
tries to seize the fundamental ingredients of the null 
world-tube of a non-expanding apparent horizon. In do- 
ing this, the world-tube is endowed with some additional 
geometrical structures that are intrinsic to the null hy- 
persurface |8| . The specific amount and nature of these 
extra structures depend on the physical problem that one 
wants to address. This introduces a hierarchy of struc- 
tures in the formalism which turns out to be very useful 
for keeping track of the hypotheses that are assumed to 
hold, as will become evident in Sec. IV. 

Before describing these structures, let us emphasize 
the change of strategy with respect to Sec. H: while 
the relevant geometry was there that of the spacelike 
three-surface where the initial data are given, the rel- 



evance corresponds now to a null three-geometry. The 
combined use of these two complementary perspectives, 
each of them suggesting their own natural geometrical 
objects, will prove to be specially fruitful. 

A. Non-Expanding Horizons 

1. Definition 

A first level in the hierarchy of structures entering the 
IH formalism is the notion of Non-Expanding Horizon 
(NEH), which incorporates the idea of quasi-equilibrium 
sketched above. We say that a hypersurface A in a vac- 
uum space-time is a NEH if 8]: 

i) It is a null hypersurface with S 2 x(lc R) topology. 
That is, there exists a null vector field on A, 
defined up to rescaling, such that g^ v l»v u = for 
all vectors w M tangent to A. The degenerate metric 
induced on A by g^ v will be denoted by q^ v . 

ii) The expansion 6m = 9 MI/ V M ^ of any null normal 
^ M vanishes on A. 3 

in) Einstein equations are satisfied on A. 

Matter can be included without problems in the scheme, 
but we will focus here on the vacuum case. 



2. Main consequences for our problem 

Let us first note that the cross-sections S ~ S 2 of the 
NEH are not necessarily strict apparent horizons since 
they are not imposed to be outermost surfaces and no 
condition is enforced on the expansion 0^ of the ingoing 
null vector (see footnote 2). Abusing of the language, 
we will however refer throughout to the cross-sections as 
apparent horizons, a practice ultimately justified in our 
problem by a sensible choice of freely specifiable data on 
the initial surface. 

a) Constant area. Owing to the null character of A, 
any null generator P defines a natural evolution on the 
hypersurface, in such a way that the area of the apparent 
horizons (a = J s d 2 V — f s yfqd 2 q, where we use the 
natural metric q a b on S induced by q^ v ) does not change, 
since £/ (ln^/g) = 6m = 0. Therefore, there is a well- 
defined notion of radius of the horizon, J?a = vW^ 71 ")- 

b) Surface gravity. Since l^ is null and normal to A, it 
can be shown to be pre-geodesic and twist free. Hence, 4 

V«i" = k {1) l» , (14) 



2 That is, a surface S in St on which the expansions 8m and 



- - (l) anu V(k) 3 v a 

of the outgoing and ingoing null vectors, P and k» respectively, Here is an y tensor satisfying q pa q p ^qau = q^u- 



satisfy Om = and 6^ < 0. 4 The symbol = denotes equality on the horizon A 
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where Km is a function on A that will be referred to as 
surface gravity (sec Appendix A). 

c) Second fundamental form 0^„ on A and evolution 
Killing vector on A. We introduce the second fundamen- 
tal form of A 16] 

= ~P a l >P f> u£iq af > = \q a ^ v dq a p , (15) 

where P a p = 6%+k a lp and q a /3 = 5^ + k a lf} + l a k Pl with 
k^ = — 1. This is an essentially two-dimensional object 
living on apparent horizons, such that 6tn = 0% while 

©ah = ^(l)Qab + o-(i) ab (16) 

defines the shear am b associated with Z M . Since Einstein 
equations hold on A, so does Raychaudhuri equation. In 
vacuo, and since the twist of l 11 cancels, it has the form 

A%) = - - (7 (i)ab a (i) ab ■ ( 1? ) 

The vanishing of 8m throughout A (so that, in particular, 

= 0) implies then the vanishing of the shear <T(;) ,. 

As a consequence, = 0. From Eq. (|T5|) we then 
see that, on A, is Lie-dragged by the null vector P. 
Therefore, although in general there is no Killing vector 
of the full space-time, the induced metric on A admits an 
intrinsic Killing symmetry. This fact extracts from the 
stronger notion of Killing horizon |17| the relevant part 
for our problem. 

d) Connection tu^. The vanishing of Mt , and the fact 
that P is normal to A suffice to define a one-form uj^ 
intrinsic to A, such that 

v'VJ* = v u lj u I^ (18) 

for any vector tangent to A. This one- form provides 
a strategy for computing in Eq. (|14fl : 

K(0 = ■ (19) 

In addition, we will see that it plays a central role in in- 
troducing the next level of the IH hierarchy of structures. 

e) Transformations under reseating of V 1 . For later 
applications, let us also summarize the transformation of 
the main geometrical objects under a rescaling of l^ by 
a function A on A. Under a change — > A we find 

q^u -> q^v , w M — > Wj, + P v In A , 

e MW AO FJ Km -» A + Z^V^A . (20) 

It is obvious from these expressions that the characteri- 
zation of NEH does not depend on the rescaling of l^. 



3. 3+1 perspective of Non-Expanding Horizons 

As discussed above, we want to cope with intrinsic evo- 
lution properties of apparent horizons. However, in con- 
trast with the previous discussion on NEH, in our Initial 
Data Problem we only dwell on a given spatial slice (at 
most, on two infinitesimally closed slices in the CTS), not 
on the whole world-tube. Therefore, we must find a pro- 
cedure to characterize an apparent horizon as a section 
of an IH by only using information on the initial spa- 
tial slice. From the NEH definition, a NEH of infinitesi- 
mal width is implemented if, together with the condition 
dm = 0, we are able to enforce = on the initial 

sphere S. The Raychaudhuri equation 117|l leads then 
to the characterization given in Ref. |1S|. that can be 
expressed as: 

The infinitesimal world-tube of an apparent horizon S 
is a NEH if and only if the shear <7(;) , of the outgoing 
null vector vanishes on S. 

Of course, if we want to extend the NEH character to a 
finite world-tube, we need to find a way to impose these 
conditions on a finite evolution interval, something that 
is not possible in the Initial Data Problem. At least, this 
instantaneous notion of equilibrium must be completed 
with a proper choice of dynamical content in the free data 
on the initial Cauchy surface. Summarizing, we see from 
Eq. (|16|l that the condition that we must impose on the 
sphere S in S t in order to have a section of a NEH is: 

Oab\s = (21) 

where the symbol | 5 stands for evaluation on S. 

B. Weakly Isolated Horizons 

A NEH describes a minimal notion of quasi-equili- 
brium, but it is not rich enough for allowing the as- 
signment of well-defined physical parameters to the black 
hole. In order to do so, we must endow the horizon with 
extra structure. Noting that the key property of the NEH 
is that l^ is a Killing vector of the metric induced on the 
horizon, a way to introduce new structure consists in en- 
forcing that other objects are Lie-dragged by Z p . 

A simple choice in this sense, that permits a Hamil- 
tonian analysis leading to (quasi-)local physical quanti- 
ties associated with the black hole, is to demand that 

Ciuj^ = 0. However, the transformation rule of w M in 
Eq. I|20l) precludes this condition to hold for every null 
normal Nonetheless, a consistent way to impose it is 
by introducing the notion of Weakly Isolated Horizon: 5 



The one-form o; M encodes some of the components of a connection 
V on A compatible with the degenerate metric lUl . This 
connection V is in fact unique as a consequence of the NEH 
definition. The stronger condition [£; , V] = defines a Strongly 
Isolated Horizon, a much more rigid structure. 
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A Weakly Isolated Horizon (WIH) is a NEH endowed 
with an equivalence class [l^] of null normals (V^ and Z** 
belong to the same class if and only if l ,fl = c Z M with c a 
positive constant) such that: 

Cilj^ = . (22) 

This condition turns out to be equivalent to (see Ap- 
pendix A) 

d(K(o)=0, (23) 

so that the zeroth law of black hole thermodynamics, 
Km = constant, characterizes the WIH notion. 

It is worth commenting that, given a NEH, it is always 
possible to select a class of null normals such that A 
becomes a WIH. Actually there exists an infinite freedom 
in the construction of the WIH structure y|. Namely, if 
the surface gravity k^ is a (non-vanishing) constant for 
a certain class of null normals [Z M ] , the same happens for 
any of the classes obtained by the non-constant rescaling 

P = [1 + 5(0, ^)e -,t «*]P, (24) 

where B(9,(f>) is an arbitrary function on S and v is a 

coordinate on A compatible with Z M , i.e., Civ = 1. In fact, 
the above rescaling does not modify the constant value of 
the surface gravity [this follows from the transformation 
rule for Km in Eq. i|2U)l ]. 

Since it is always possible to find WIH structures on 
a given NEH, the WIH concept does not correspond to 
a real restriction on the physics of the system. However 
it does impose a restriction on the spacetime slicing by 
the hypersurfaces St introduced in Sec. [H] if we tie l^ 
to t, i.e. if we impose that there is a member f of the 

WIH class such that Lit = 1. We call such a slicing 
a WIH- compatible slicing. 

The derivation of the mass and angular momentum 
expressions for a WIH using Hamiltonian techniques is 
beyond the scope of this work (see Refs. 00). Here, 
we will simply extract those points which are relevant for 
our analysis. The general idea is to characterize physical 
parameters as conserved quantities of certain transforma- 
tions that are associated with symmetries of the WIH. A 
vector field tangent to A is said to be a symmetry 
of the particular WIH under consideration if it preserves 
its equivalence class of null normals, the metric q^ v , and 
the one-form , namely 

CyF = constant • Cyq^v = , Cycu^ = . (25) 

In Sec. IV we will be interested in non-extremal black 
holes, for which Km ^ 0. In that case, the general form 
of a WIH symmetry is 

= cvl" + b v S^ , (26) 

where cy and by are constant on A and is an isometry 
of the apparent horizon S. 



The definition of the conserved quantities goes first 
through the construction of an appropriate phase space 
for the problem and then through the analysis of canoni- 
cal transformations on this phase space ■ An important 
point is that the relevant transformations are generated 
by diffeomorphisms in space-time whose restriction to the 
horizon A are symmetries of the WIH in the sense of l|25fl . 

1. Angular momentum 

In order to define a conserved quantity that we can 
associate with a (quasi-)local angular momentum, we as- 
sume that there exists an azimuthal symmetry on the 
horizon A (actuall y, t his hypothesis can be relaxed; see 
in this sense Ref. [l5|). Therefore, we assume the exis- 
tence of a vector ^ tangent to S C A, which is a 5*0(2) 
isometry of the induced metric q a b with 2tt affine length. 

The conserved quantity associated with an extension 
of ipV to the space-time is given by [7| 

Ja = -^J s (Av d 2 V = ^ J/^ K ^ d2y > ( 27 ) 

where for convenience we have expressed it in terms of 
objects in the 3+1 decomposition. In particular, s l is the 
outward (pointing towards spatial infinity) unit vector 
field in S t normal to the apparent horizon S. 

2. Mass and boundary condition for t^ 

The definition of the mass is related to the choice of 
an evolution vector t^ with appropriate boundary con- 
ditions, namely that t^ — > (dtY L at spatial infinity and 

_> - Q t¥? M with P € and Q t constant on the 
horizon (note that t^A is a WIH symmetry). The deter- 
mination of the mass expression proceeds in two steps. 

First, the vector has to satisfy certain conditions 
to induce a canonical transformation on the phase space. 
This turns out to be equivalent to the first law of black 
hole thermodynamics 0, , whose practical consequence 
for us is that the mass M, the surface gravity Km, and 
the angular velocity f2 t depend only on the radius Ra 
and the angular momentum J a of the black hole, 

M = M(R a ,Ja) , 

fit = ntOR A ,JA) , (28) 

but without determining the specific functional form. It 
is worth emphasizing that this dependence on Ra and 
J A (though arbitrary in principle) must be the same for 
all solutions to the Einstein equations containing a WIH. 

In a second step, this dependence is fixed to coincide 
with that found in the stationary Kerr family of black 
holes. This is not an arbitrary choice, but a normaliza- 
tion consistent with the stationary solutions. Technically, 



6 



this is accomplished by requiring that, at the horizon A, 
t^ + Vt t tp^ reproduces just the null normal (in the consid- 
ered class [I 11 ]) whose surface gravity equals that of the 
Kerr case, something that is always possible in the non- 
extremal situation via a constant rescaling. This singles 
out a vector t% , satisfying 



with 



ta+Q Kerr (R A ,j A )^^ci" = ^ 



K K „ T (Ra, Ja) 
K(I) 



(29) 



(30) 



as the evolution vector used for the derivation of the mass 
formula. The final expressions obtained in this way for 
the physical parameters of the horizon are: 



M A = 

«A = K Kerr (RA,JA) 

A = & Kerr (RA, Ja) 



AG 2 Jl 



2GR A 
R A — 4G 2 J A 



2R 3 Ay /R A + AGJl 
2GJ A 



R A yjR\ + AGJ'i 



(31) 



IV. DERIVATION OF THE BOUNDARY 
CONDITIONS 

We are now in an adequate situation to derive bound- 
ary conditions for the elliptic equations in Sec. II. In 
doing so, we adopt a coordinate system {t^x 1 ) stationary 
with respect to the horizon, in the sense that the null tube 
A can be identified as the hypersurface r(x l ) = constant 
for a certain function r which is independent of t. It can 
be shown that this happens if and only if i M is chosen 

tangent to A, i.e. l^ = 0. 

Since we want to have a notion of angular momen- 
tum for the black hole, following the discussion in Sub- 
sec. III.B.l we make the hypothesis that our physical 
regime permits the imposition of an axial isometry tp a on 
S ~ S 2 C A. Even though this is a strong physical hy- 
pothesis (especially when having in mind binary black 
holes), we must emphasize that the bulk space-time will 
still be generally dynamical in an arbitrarily close neigh- 
borhood of the horizon and that (p^ does not need to 
extend to an isometry there. 

To construct the equilibrium black hole on S, we follow 
the steps dictated by the hierarchy of the IH formalism. 



A. Adapting the evolution vector to the horizon 



to the horizon by relaxing to a NEH the particular WIH 
structure implicit in Eq. J2HJ)- That is, we only impose 



t" 



n A ^ ~ j" 



(32) 



where the proportionality needs not be given by a con- 
stant on A. Using the proportionality 



(33) 



(where s M is again the outward unit spatial vector normal 
to S) and the decomposition (JJJ of i M in terms of the lapse 
and shift, we conclude 



/3* = as* - n A tp* 



(34) 



from which boundary conditions for the shift on S im- 
mediately follow. 

Actually the choice of stationary coordinates with re- 
spect to the horizon automatically leads to the expression 

f3 l = as 1 — W l , where W l is an arbitrary vector tangent 
to S. So, here we enforce W l to be precisely Sl A tp l , a 
choice that will in fact simplify the imposition of the 
NEH structure. 6 



B. Non- Expanding Horizon condition 

We now properly impose the NEH condition. As men- 
tioned in Subsec. III. A. 3, in this Initial Data Problem 
we demand S to be a slice of a NEH of infinitesimal 
width. For this, we impose condition 121|) . Owing to the 
rescaling property (|20(l of MI , under an arbitrary (not 
necessarily constant) rescaling of / M , and taking advan- 
tage of the t^ adaptation to the horizon implemented by 
the shift boundary conditions, we can write 



(35) 



In our stationary coordinates with respect to the horizon, 
this simply reads 



= 2Q ab = d t q a b + C v q ab 



(36) 



But, under our hypothesis about the existence of an axial 
isometry on the horizon, the second term must vanish on 

its own: C v q a b = 0. Using q ab = j ab (for angular co- 
variant components) we find 



dtlab = 



£ v 7afc = 



(37) 



In particular, the restrictions must hold on S. These are 
the NEH boundary conditions. Note that their simple 



Aiming at imposing the NEH structure, but already 
motivated by the boundary condition for the evolution 
vector selected by the determination of t%, we adapt i M 



6 Although the derivation of Ja actually involves a WIH structure, 
its expression can be shown to be already well-defined for a NEH. 
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form depends critically on the specific choice made for 
W % in the previous subsection. 

Using now the conformal decomposition of the metric, 
these conditions translate into 



0, 



(4% b £ v y + ^£ v j ab )\ s = Q. 



(38) 
(39) 

The crucial feature, and the ultimate reason for using the 
CTS, is that these conditions can be satisfied by an ap- 
propriate choice of the free data ^ a b and dt"/ a b- Condition 
H39J) . expressing the axial symmetry of the horizon, must 
be enforced by a self-consistent selection of the free data 
7ah on S (^f is a functional of 7y). Regarding (J2HJ, if we 
first take its trace with respect to the conformal coun- 
terpart of the metric q ab induced on S, q ab = ^~ i q a b 
(satisfying q acx i c b = 5%) and then use Eq. ijTU|) for 9 t 4", 



A -„ 



we find (calling tr s j = q ab d t q a b = q dtlab) 



g" 



aK + -tr s 7 







(40) 



In addition, from Eq. I|38[l it follows that the g-traceless 
part of <9t7ab must vanish. Therefore, on the boundary 
S, this part of the free data has the form 



dtlab - 2^57 lab 



= 0. 



(41) 



Condition l|4U|) is an inner boundary condition for . 
Since we have imposed on S the Dirichlct boundary con- 
ditions Ij34(l on j3 l , we have no direct control on the sign 
of Di(3 l there. In order to guarantee the positivity of W 
via the application of a maximum principle, the factor 
multiplying W in Eq. I|40|l should be non-negative. The 
analytical study of this issue goes beyond the present ge- 
ometrical derivation (see Ref. 0] for a discussion on this 
point in a related context). We simply comment that the 
choice of free data for tr s 7 [and, more indirectly, that of 
the radial components of the free data 7 y , which deter- 
mine s l — 7 r VV7 w m condition l|34|) ] could play a key 
role in ensuring that Eq. 1)41) [) is a well-posed condition. 

Except for the implicit use of a well-defined concept of 
angular momentum based on the WIH formalism, the no- 
tion of quasi-equilibrium provided by the NEH structure 
has proved to be sufficient to set boundary conditions for 
the Initial Data Problem, since it prescribes boundary 
values for (Hamiltonian constraint) and (3 l (momen- 
tum constraint). If this is the problem that we want to 
solve, we can stop here. However, if we want to solve 
also the trace of the evolution equations, we need to find 
appropriate boundary conditions for the lapse. We will 
show how the existence of a WIH structure can be ex- 
ploited with that aim. 



isolated. In fact, the determination of this class is not 
unique, but there exists an infinite freedom of choice. 
In this subsection, we will first discuss the restrictions 
on the lapse function that follow from the introduction 
of a WIH-compatible slicing and then employ the free- 
dom in the choice of WIH to suggest possible boundary 
conditions for the lapse that are specially suitable for 
numerical integration. 

Let us start by choosing 1% as the representative of the 
class of null normals for the WIH, [l^]. The inner bound- 
ary condition 129|) employed in the determination of the 
mass formula then singles out an evolution vector t% on 
the horizon. We proceed as in Subsec. IV. A, but impos- 
ing to coincide exactly with therefore demanding 
that the (E t ) foliation constitutes a WIH-compatible slic- 
ing. According to the characterization 12: ill of the WIH 
notion, the surface gravity is constant. We further 
assume ^ 0, thus restricting the analysis to the non- 
extremal case. We write 1% in terms of 3 + 1 objects, 



l» = aP = a{n^ + s") 



(42) 



where we have explicitly defined the vector P. Again, 
the introduction of the lapse and shift decomposition for 
t% in Eq. i|29|) leads to the boundary conditions l|34|) for 
the shift. In order to analyze the conditions on a, we 
calculate the expression for Kn\. We proceed in several 
steps: 

1. Contracting Eq. (|18|) . particularized to the one- 
form u p associated with Z^, with the ingoing null 
covector k p = (n^ — s M )/2 and expanding the result- 
ing expression, we find for any vector tangent to 
A, 



(43) 



Employing the definition of the extrinsic curvature 
(0 and the identity n p V p n v — r ) p u V p In a, we get 

w"w„ = -i/V [K^ + rvtfV, In a) . (44) 

2. Taking I 1 * as the tangent vector v p and remember- 
ing expression (|19f) . we obtain 

= f"w M = s"V M In a - s»s v K^ . (45) 

3. Recalling then the transformation (|2U|I of the sur- 
face gravity under a rescaling of the null normal, 

k {1o) ^a(s p V p Ina H-i^lna. (46) 



C. Weakly Isolated Horizon condition 

As we have commented, given a NEH one can always 
find a class of null normals so that it becomes weakly 



4. Finally, imposing that Ka ) equals K Kerr (RA, J/\), 
we find 

h k „,XRa, Ja) = s i D i a - s i sPK i j a + d o In a . (47) 
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This restriction, arising from the WIH-compatible slic- 
ing condition, can be regarded as an evolution equation 
for the lapse on the horizon. Properly speaking, it is 
not a boundary condition for the initial data, because it 
contains the derivative of the lapse in the direction of 1%. 

Actually, by exploiting the freedom of choice in the 
WIH structure, one can freely set the value of the lapse on 
the initial section of the horizon, iS. This is a consequence 
of the fact that, from Eq. I|24|) and the relation 1% — od^, 
a change of WIH structure results in a rescaling of the 
lapse: 



1 + B(6>,^)e~ K ^"'( iW A)t 



(48) 



Since S can be identified, e.g., with the section t = of 
A, the initial value of the lapse on the horizon gets mul- 
tiplied by an arbitrary positive function on the sphere. 
Therefore, it can be chosen at convenience, at least as 
far as the demand for a WIH-compatible slicing is con- 
cerned. In fact, if we were solving an evolution problem 
by following a constrained scheme (see e.g. Ref. 
this initial choice for a together with Eq. (|47|) might be 
used to set inner boundary conditions for the lapse at 
each time step. 

Once realized the freedom in the choice of the initial 
value for a on S that follows from the dynamical charac- 
ter of Eq. (|47|l on A, one may ask whether it is possible 
to benefit from this arbitrariness and put forward a par- 
ticular proposal for the choice that could be considered 
specially advantageous. In this sense, one would like to 
ensure that, under evolution on A, the lapse will neither 
increase exponentially nor decrease to (zero or) negative 
values. Apparently, the best way to favor this, at least 
locally, is to pick up, among the infinite WIH structures, 
that in which the Lie derivative of the lapse with respect 
to the null normal vanishes initially: (Ci o a)\ s = 0. Using 
Eq. 148|) . one can prove under very mild assumptions 
that such a choice of WIH structure exists. Adopting it, 
Eq. (|47|l becomes a true boundary condition for the lapse 
on S: 



[s l D l a - s*s j Kij a] | = K Kerr (i? A , J a) ■ 



(49) 



Note in fact that, to deduce this condition, one only 
needs to demand [in the passage from Eq. to Eq. 

(P7J) ] that K( /o j coincides with the constant n Kcrr (Ra, Ja) 
on S, and not on the whole of A, because one finally re- 
stricts his attention just to the initial section of the hori- 
zon. As a consequence, and in contrast with the situation 
found for the NEH conditions, the above prescription for 
the lapse on the boundary is only a necessary condition 
for specifying a WIH of infinitesimal width. The extra 
condition that be constant in the rest of A, namely 
(£/ o K(; o ))| s — 0, would involve the evolution equations 
and the second time derivative of the lapse, and there- 
fore cannot be imposed in terms of the initial data. 

Finally, we comment that an alternative way of dealing 



with the WIH condition k 



constant would consist in 



interpreting Eq. (|47() as a constraint on the free data on 
the inner boundary. 



D. Binary quasi-circular orbits 

In the previous subsections we have characterized the 
quasi-equilibrium state of each horizon exclusively in lo- 
cal terms. However, the study of a binary black hole in 
quasi-circular orbits requires, in addition, a global notion 
of quasi-equilibrium. In the general case, such a global 
quasi-stationary situation is described by the existence 
of a global quasi-Killing vector L^. In the binary black 
hole case, this is a helical vector (see Refs. |l(J,|2(|) that 
Lie drags the horizons, i.e. L^Ia is tangent to each hori- 
zon A. Imposing asymptotic flatness, we have at spatial 
infinity L M — > t(£, + f2 or (, , where t{£, and are vec- 
tors associated with an asymptotic inertial observer and 
Qorb is the orbital angular velocity. We can adapt the 
coordinate system, introducing an evolution parameter 
t' such that I/ M = (<9f) M . In such a case, with the 3+1 
decomposition — a'n^ + outer corotating bound- 
ary conditions follow 

lim P H w n orb 0^ , lira a' w 1 , lim*wl. (50) 

r — >oo r — >oo r — >oo 

In these coordinates, one chooses the time derivative part 
of the CTS free data to vanish 



d t >% = , d v K = 



(51) 



In the general case, L m |a and t% defined in Eq. (|29|) do 
not coincide. Since they are both tangent to the horizon, 



^ = p K 



x' 



(52) 



where p is a scaling factor and x M is tangent to A with 
= 0. As a consequence, if we adapt coordinates 
to L^, hence using t' , the expressions given in the previ- 
ous subsections must be corrected. We will discuss two 
possibilities. 

a) Corotating coordinate system (fully adapted to L^). 
From Eqs. I|29|l and (|52|l we can write 



L^ + p n A ip" - X 



" = p 



(53) 



A natural ansatz for is given by fl or b where M 
is the azimuthal vector tangent to each horizon and as- 
sociated with the normal direction to the orbital plane. 
Since Qa provides a well-defined notion of rotation angu- 
lar velocity, we can define the corotating physical regime 
in an intrinsic way as the case with p = 1 and Q or b = Ha, 

from which = l£ follows. More generally, proceeding 
as in Subsec. IV. A we find 



[3 H = a's 1 



pflA 



X 



(54) 



Imposing the axial symmetry on each horizon A, we 
deduce again condition (|39|l . Defining 



choosing a priori the values of a and Ci o a on S and then 



"if 



p flA ^ 



(55) 
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the requirements 6 a h = and dfjij 
lead then to the conditions: 



[see Eq. JST}] 



1 

8*3 • 



0. 



= 0,(56) 
(57) 



which replace Eqs. I|4(J|) and (|41|l . respectively. Finally, 
the condition for the lapse is still derived as in Subsec. 

IV. C. From Eq. it follows that a' = pa, and hence 
it can be shown 



sWia' - s i s j K ij a' 



K Kerr (i?A, Ja) 



(58) 



Convenient ansatze for p and x^ 1 must be introduced in 
practice to cope with these conditions. 

b) Warped coordinate system. An alternative choice 
consists in adopting an evolution vector t^ such that its 
boundary value on each horizon coincides with t%, but 
adapts itself to L M at a typical distance S from them (in 
this sense, the distance between the black holes provides 
a natural length scale in the binary problem). Hence, 
this vector t^ interpolates between t% and L M , warping 
the coordinate system to better accommodate the phys- 
ical situation in each of the considered spatial regions 
(note that the vector L M follows the translational mo- 
tion, whereas t% is adapted to the intrinsic rotation on 
the horizon). Of course such a coordinate system can re- 
main regular only for a finite amount of time (typically 
one orbital period). 

In practical terms, this coordinate system is defined 
by the outer boundary conditions J5UJ, without primes 
in (3 l and a, and the inner boundary conditions (|34() [on 
S] , (|40|l , (|41|l , and (|49|l on the constrained and free data. 
Moreover, the function ir s 7 does not have to vanish on 
S, therefore helping to ensure the positivity of *, even 
though it must be negligible at a distance of order 5. 
Likewise, dtK and the radial components dtj r j become 
roughly zero at a distance d of each horizon. Thus, in 
this coordinate system, the conditions on the horizons 
are easier to impose, there is no need to worry about the 
factor p, and one gains control over the positivity of *. 

Once the time derivative part of the free data has been 
fixed, either in the corotating or in the warped coordi- 
nate system, one would have to consider the rest of the 
free data. The choice of the conformal metric must be 
consistent with restriction (|39[1 [and l|57|l in corotating 
coordinates] and subject to the constraint det( 7 jj) = /. 
The adequate determination of the physical content of 7y 
goes beyond the limited scope of this paper and must be 
addressed by means of a proper analysis of the stationary 
regime for the evolution equations J5j- 

As for gauge fixing, the Dirac gauge in Ref. Q ap- 
pears to be a quite natural choice for the spatial one in 
the CTS setting, whereas the boundary condition 14UI) 



might be viewed to suggest a maximal slicing (K = 0) 
for the temporal gauge in order to improve the control on 
the positivity of *. However, this latter gauge is not com- 
patible with coordinates of Painleve-Gullstrand or Kerr- 
Schild type, which are actually appropriate for the shift 
boundary condition l|34|) . We do not here subscribe to 
a particular fixation of the gauge, allowing an optimal 
adaptation to each case considered. 



V. COMMENTS ON PREVIOUS APPROACHES 



A. Cook's 2002 proposal 

Inner boundary conditions for the elliptic Eqs. Qllfl- 
l|12p. and in the quasi-circular regime of a binary 
black hole system were presented in Ref. . The scheme 
proposed in that work starts by imposing on each excised 
sphere S the presence of a Killing horizon, together with 

an apparent horizon condition, 9>j\ = 0, where f M is de- 
fined by Eq. jlSjl. Denoting by C,^ = an^ + j3 v ' s v s^ the 
component of i M orthogonal to S, the following quasi- 
equilibrium conditions were imposed: 

1. The inner boundary S remains an apparent hori- 
zon: £>C,0(h = 0. 

2. The expansion 9^ associated with the ingoing null 
vector A; M = (rt M — s M )/2 does not change in time: 
W(k) = • 

These conditions are enforced under the approximation, 
motivated by the stationary case, that the shear <7/a as- 
sociated with the outgoing null vector vanishes. 

Before we actually compare the resulting boundary 
conditions (or rather some ellaboration of them; see the 
last part of this section) with those of Sec. IV, we make 
some general remarks on the involved quasi-equilibrium 
conditions. 

Under the vanishing shear approximation, the condi- 
tion C^9,i[\ = leads to (i v s v = a, thus making ^ a 

null vector parallel to Z M . In particular, this implies that 
the underlying coordinate system is stationary with re- 
spect to the horizon. Therefore, this condition is either 
redundant with the vanishing shear approximation (via 
Raychaudhuri equation) or must be considered as a gauge 
choice, and not as an actual quasi-equilibrium condition. 

More generally, in Ref. the conceptual status of the 
vanishing shear hypothesis is not clearly stated and an 
explicit prescription for imposing it in terms of the initial 
data, such as Eqs. (|41|l or (|57|l . is missing. The IH analy- 
sis shows that the vanishing of the shear is the key quasi- 
equilibrium condition: it guarantees that the world-tube 
of apparent horizons is a null hypersurface. More ex- 
plicitly, if o~ /j\ = is not taken as a quasi-equilibrium 
characterization but only as an approximation that might 
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occasionally fail, the vector is no longer necessarily 
null. As a consequence £>c,9rt\ = would not really be a 
quasi-equilibrium condition (for instance, C^6,j\ vanishes 
also for dynamical horizons |2lJ, where £ M is spacelike). 

Hence, as already pointed out in Refs. [THlll3 |. the ap- 
proach of Ref. [9| is very close in spirit to that encoded 
in the IH formalism; in fact, if the approximation of van- 
ishing shear is eventually satisfied, a NEH is actually 
constructed. However, its quasi-equilibrium conditions 
can be refined 7 (see Ref. [23 and Subsec. V.A). By 
contrast, a virtue of our approach, fully based on the IH 
scheme, is a clear identification and understanding of the 
physical and mathematical hypotheses that characterize 
the horizon quasi-equilibrium. 

The need to clarify, from a conceptual point of view, 
the quasi-equilibrium hypotheses in Ref. |0j can be illus- 
trated as follows. In Ref. j^, the boundary conditions 
derived in Ref. 9], identified as IH conditions, are dis- 
regarded as technically too complicated. They are then 
substituted by a heuristic set of conditions, involving in 
particular = c^lny^ = Dif3 z — aK. This condition, 
which is equivalent to = 0, turns out to be the 
NEH condition on the conformal factor in the corotat- 
ing physical regime [make 77^ = in Eq. i|56|) ]. under the 
quasi-equilibrium bulk condition d t jij = assumed in 
Ref. |23. Besides, the Kerr-Schild data (motivating the 
boundary values for the shift in that reference) are con- 
sistent with Eq. (|54(l . At the end of the day, we find that 
the heuristic choice turns out to be one which is truly in 
the spirit of the IH scheme. Let us nonetheless mention 
that these boundary conditions are not imposed on the 
horizon itself, but in its interior. 



B. Addendum 

After the first submission of this work, a paper by Cook 
and Pfeiffer appeared [24[ which provides a refinement of 
the discussion and proposals made by Cook in Ref. @. 
We now comment on the relation between our approach 
and the quasi-equilibrium and boundary conditions pro- 
posed in Ref. |2J] (and in Q), in order to facilitate the 
comparison of our results with those of that reference. 

1. Quasi-equilibrium conditions. Quasi-equilibrium is 
characterized in Ref. [24| by the geometrical conditions 



'(0 



= 0, which are exactly those required to 
construct a NEH horizon, as discussed in Subsec. IV. B. 
Since no other condition is imposed (the requirement 

^C^(fc) = of Ref. 9] is dropped), the analysis remains 
at the level of a NEH in the IH hierarchy, whereas our 
approach explores the WIH structure. 



2. Condition on . In Refs. H 13, this boundary 
condition follows from the requirement of vanishing ex- 
pansion for an apparent horizon. It is therefore es- 
sentially equivalent to Eqs. I|40|) and l|56|) . However, the 
mathematical expression derived from 9^ = adopts 
different forms [see also Eq. <|A13|) in Appendix A]. 

3. Condition on [3 l . In corotating coordinates, the 
boundary condition (79) of Ref. 9] essentially coincides 
with our Eq. (|54() . A crucial refinement is introduced in 
Ref. [24| by actually imposing that the shear vanish: the 
projection of the shift on S [our vector —rf in Eq. I|55|l ] 
must be a conformal symmetry of q a \,. This is equivalent 
to our condition ijSTjl (see also Appendix A). 

The main difference between both approaches is our 
demand of an azimuthal symmetry ip^ for the metric q a b, 
namely Eq. i|39H . On the one hand, this makes condi- 
tions in Ref. [24| more general than ours but, on the 
other hand, thanks to this symmetry we are able to in- 
troduce a definite, intrinsic spinning angular velocity f^A 
which, together with £l or b, permits to analyze the rota- 
tional regime of the system (corotational, irrotational or 
general case). 

In addition, the availability of naturally leads us to 
consider Eq. (|54|l as a boundary condition on the shift. 
As a consequence, Eq. I|57f) becomes a constraint on the 
free data "f a b and dt"f a b, rather than providing a boundary 
condition for rj^ as in Ref. |24j . 

4. Condition on a. The analysis of a WIH carried out 
in Subsec. IV. C shows that the initial boundary value for 
the lapse is basically free. This conclusion is also reached 
in Ref. [24| after a numerical study. It is worth discussing 
the relation between the proposals that have been made 

for the choice of this boundary value. Condition C^9^ = 
in Ref. 9] can be written as 



Va 



(59) 



with V defined in Eq. (85) of Ref. Our requirement 
(|l9*|l and Eq. (|59*|) are simultaneously satisfied only if 
on S. This is a nontrivial identity, so 



Va 



7 We have obviated the discussion of Ref. |S| about the way to 
enforce the horizon to remain in the same coordinate location, a 
discussion that can also be simplified. 



that both conditions are generally different. 

Insight on their relation is provided by Ref. Q , where 
the freedom in the construction of a WIH structure [Z M ] 

is fixed by imposing that Cn\Qn.\ = (with k^l^ = — 1), 
once it is assumed that a certain operator M which acts 
on S has a trivial kernel [see Eq. (4.8) of Ref. @ for 
the definition of M and note its close connection with 
V}. This analysis can be applied to study the possible 
degeneracy of condition l|59|l in terms of the invertibility 
of M if, in addition, it is satisfied that k,r is constant. 
If that is the case, employing Eq. I|45(l one can check 
that conditions 1)49(1 and l)59|l coincide only if the lapse 
is constant on the boundary. More details on this issue 
will appear elsewhere. 

Notice that the choice of representative made in a 
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WIH class via the lapse boundary condition i(49() deter- 
mines the initial lapse once 7.^ and K+j are given. In 
this sense, the condition for the lapse is not problem- 
atic by itself. However, our full set of boundary condi- 
tions for a, iff and f3 l , together with the choice of free 
data (jij , dtj' 1 -' , K,dtK), may not be sufficient to single 
out a unique solution to the Initial Data Problem. In 
fact, this degeneracy seems to occur when our boundary 
conditions are implemented in the spherically symmetric, 
time-independent case if one uses a maximal slicing and a 
flat conformal metric 8 . Nonetheless, the presence of this 
degeneracy may depend on the actual choice of initial 
free data (e.g., isotropic coordinates in the commented 
example). For each specific choice, it is generally only 
after a numerical study that one may decide whether a 
degeneracy exists. 

VI. CONCLUSIONS 

In this work we have explicitly shown how the IH for- 
malism provides a rationale for some aspects of the nu- 
merical construction of initial data for a space-time con- 
taining a black hole in local quasi-equilibrium, with spe- 
cial emphasis in the binary case. 

The IH framework sheds light into the justification and 
implications of already existing quasi-equilibrium sets of 
conditions for the analysis of this problem. The hierar- 
chical structure of the IH formalism permits a control on 
the hypotheses that arise at each of the considered steps. 

Adopting the IH approach fully, we have derived a set 
of boundary conditions on each black hole horizon for 
solving the elliptic equations obtained in a CTS scheme, 
deduced from the constraints and the trace of the evolu- 
tion part in Einstein equations. 

In a first step, the NEH condition characterizing quasi- 
equilibrium (O a fc = By) = cr(i) ab = 0), together 
with the choice of spatial coordinates stationary with re- 
spect to the horizon, provides boundary conditions for 
the shift [see Eqs. 1(34(1 and 1)54(1] and the conformal fac- 
tor [see Eqs. (|40f) and l|56|)]. These conditions are ba- 
sically equivalent to those of Ref. |9j (at least in the 
recently refined form presented in Ref. |2J|). In a sec- 
ond step, the requirement for a WIH-compatible slicing 

(«(i) = constant ^ 0) leads to the evolution equation i|4"7|) 
for the lapse on A, leaving the choice of its initial value 
essentially free. Once this point has been acknowledged, 
we have tentatively suggested a specific boundary con- 
dition for the lapse in Eq. (|49|l . obtained by fixing the 
freedom which is available in the construction of a WIH 
structure. In addition to these boundary conditions, the 
NEH requirement entails that the free data on the initial 
slice fulfill, on the horizon, the constraints (|41|l [or 1(57(1 



We thank G.B Cook for pointing out this fact that also happens 
with the boundary conditions of Ref. [jj . 



in corotating coordinates] and 139|) (assuming an axially 
symmetric horizon; see Ref. |25| otherwise). Similarly, 
Eq. 1)47(1 could alternatively be seen as a constraint on 
the free data if one decided to fix the lapse on A. 

These inner boundary conditions and constraints are 
sufficient conditions for constructing a black hole in in- 
stantaneous quasi-equilibrium. However, in order to ob- 
tain black holes in quasi-equilibrium during a finite evo- 
lution time rather than just instantaneously (as required 
in the quasi-circular binary black hole problem), these 
conditions must be complemented with appropriate free 
data that encode the desired dynamical behavior. 
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APPENDIX A: SOME TECHNICAL DETAILS 

In this appendix we explain some calculations and for- 
mulas employed in the main text. 

a) Proof of Eq. ((14(1 . If A is defined as the hypersur- 
face r = constant, its normal takes the form — AV^r 
for certain function A. Then V[ M Z„] = A /r^V^i A. Con- 
tracting with Z M , we find 

I"V M Z„ = l^^hx\l v . (Al) 

Eq. (fMfl follows by identifying = V M In A. 

b) Proof of Eq. ((18(1 . Choosing the normalization of 
the outgoing and ingoing null vectors and fc M on A so 
that k^l^ — —1, the induced (degenerate) metric on A 
can be written as 

c lnv = 9iw + k)J>v + l^kp = g^u + n^riu — s^s„ . (A2) 

Expressing Ciq^ in terms of the connection V M and using 
Eq. I(A2() . we get for the second fundamental form 115() 
the formula: 

©a/3 = (fctffF^v ■ (A3) 

Expanding now q^ ' a q" ' g and recalling Eq. I(14|) . we find 

V Q //3 = 6 Q/ 3 — l a k u V v li3 

- (IfValv + lJPWVulJlfi . (A4) 

Since, on a NEH, one has Q a g = (see the text), con- 
traction with a vector w M tangent to A (so that v^l^ = 0) 
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leads to Eq. (|18fl after defining 

L0 a — (k V a_lv ~t~ l a k^k V /ill* ) 



-P^ a k»V^ . (A5) 



c) Proof of Eq. (|23[l . To demonstrate this equation 
starting from condition (|22|l . the key remark is the pro- 
portionally between the exterior derivative of w M and the 
volume 2-form on the sphere S 2 ~ S, 



duj oc y/q d 2 q 



(A6) 



This is a nontrivial result that follows from the definition 
of NEH, and we refer the reader to Ref. [|| for its proof. 
As a consequence, du lives on S and its contraction with 

the null normal of A vanishes: l^W^ui^ = 0. Therefore, 
using the Cartan identity, 

= C t u v = rv^j + V M {Vu) v ) = V,^,) . (A7) 

d) General expression of Q ab . In order to enforce the 

condition Q ab = 0, we have made use in the text of 
a coordinate system that is stationary with respect to 
the horizon, so that = + with W M tangent 
to S (W 11 is the black hole surface velocity introduced 
by Damour |16l l2q1. Then, we have considered several 
specific choices for on A, namely, the vector tJJ sin- 
gled out in Eq. (|29|) for the warped coordinate system, 
or the quasi-Killing vector for corotating coordinates. 
For the sake of completeness, we now provide the gen- 
eral expression of the second fundamental form on A for 
arbitrary vectors and VF M . 

From the definition of M1/ and Eq. (| A2|) we find 

1 



9 Qb = | (£ t q a b + 2 D a W b + 2 D b W a ) 

2 (^tQab + Z-WQab) ■ 



(A8) 



Here, 2 D a denotes the connection compatible with the 
metric q a b induced on the sphere S and W a = q a bW b . 
Note that Eqs. and %$jc$ follow straightforwardly 



from this when one imposes VT M to be the g Q fc-isometry 
The conformal decomposition q ab — ^> 4 q ab leads to 



e afc 



* 4 



1(i)q~ab + £tq~ab - t; (A Ing) q a b 



2 D a W b + 2 D b W a - 2 D C W C q al 



(A9) 



and 



e {l) = -£ t \nq + 4£ t \n^+ 2 D a W a +4W a 2 D a lnt>, (A10) 

where 2 D a is the connection compatible with q ab and 
W a = q a bW h . In particular, with the notation rf for 
and choosing i M to be the quasi-killing vector (so 
that d t q ab = 0), the vanishing of the expansion in Eq. 
QAI 0(1 leads to condition [after substituting Eq. (fTOf) 
for dt^}. On the other hand, from Eq. (|A9|) . the trace- 
less part of Q ab = is equivalent to condition (|57|l . It is 
then clear from Eq. (|A9|) that in corotating coordinates 
the vector rf must be a conformal symmetry generator 
of q a b- 

Finally, a more general expression for Qj a can be ob- 
tained if we do not assume a coordinate system stationary 
with respect to A. Substituting l^ — a(n^ + s^) in Eq. 
I(A3(I and expanding the derivative we find 

Q a 0=a(D„s v -K„ u )q^ a q u . (All) 

Taking the trace, the standard expression for 6ns follows, 

6 m = a(As i + KysV - K) . (A12) 

Hence, after a conformal decomposition, 



W = [As'D, In * + + ^- 2 KijtP - f 2 K 
xar 2 , 

where s l — ^> 2 s l . 



(A13) 
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